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Abstract 

A new mechanism for nonmodal amplification of stochastic disturbances in elasticity-dominated 
channel flows of Oldroyd-B fluids is reported. For streamwise-constant flows with high elasticity 
numbers fi and finite Reynolds numbers Re, we establish the 0{fiRe^) scaling of the steady-state 
variance amplification. This demonstrates that large velocity variances can be sustained in weakly 
inertial/strongly elastic stochastically driven channel flows. The mechanism examined here provides 
a possible route for a bypass transition to 'elastic turbulence' and might be exploited to enhance 
mixing in microfluidic devices. 

PACS numbers: 47.20.Gv, 47.20.Pc, 47.27.Cn 
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Fluid motion is usually classified as either laminar or turbulent; flows that are smooth 
and ordered (laminar) may become complex and disordered (turbulent) as the flow strength 
increases. This process is commonly referred to as transition (to turbulence). A laminar 
flow is characterized by smaller viscous shear stresses and velocity gradients than a turbulent 
flow, leading to a lower skin-friction drag force on a vehicle moving through a fluid in the 
laminar regime. On the other hand, a turbulent flow provides excellent mixing which is, for 
example, required in chemical reactors, heat exchangers, and combustion engines. 

The classical approach to transition examines the linearized equations for exponentially 
growing normal modes. The existence of these unstable modes implies exponential growth of 
inflnitesimal perturbations to the laminar flow, and the corresponding eigenfunctions identify 
flow patterns that are expected to dominate early stages of transition. This approach agrees 
with experiments in many relevant flows (e.g., those driven by thermal and centrifugal 
forces [1]) but it comes up short in matching experimental observations in wall-bounded 
shear flows (flows in channels, pipes, and boundary layers). Failure of classical hydrodynamic 
stability analysis is attributed to the nonnormal nature of the linearized equations which 
may manifest itself by transient growth of perturbations |2J, protrusion of pseudospectra 
to the unstable regions [Ij, and large receptivity to ambient disturbances |3], Hj. Even in 
stable regimes-owing to nonnormality-perturbations that grow transiently before decaying 
due to viscosity can be conflgured, irregularities in laboratory design can lead to instability, 
and background disturbances (such as free-stream turbulence or surface imperfections) can 
be amplifled by orders of magnitude. These conclusions can be reached by performing 
transient growth, pseudospectra, or variance ampliflcation analyses |5l Ej. All of these 
methods demonstrate the importance of streamwise-elongated flow patterns of high and low 
streamwise velocity in transitional wall-bounded shear flows of Newtonian fluids; this is at 
odds with classical modal stability results, but in agreement with experiments conducted in 
noisy environments [7]. 

Transition to turbulence in viscoelastic fluids is important from both fundamental and 
technological perspectives [8J. The observation that transition can occur even when the 
effects of fluid elasticity dominate those of fluid inertia-which is a primary cause of tran- 
sition in Newtonian fluids-is particularly intriguing. Improved understanding of transition 
mechanisms in viscoelastic fluids has broad applications in physics, ranging from deeper 
insight into order- disorder transitions in spatially extended nonlinear dynamical systems to 



2 



enhanced mixing in microfluidic devices through the addition of polymers [9l [TOl [TT] . 

Amphfication of stochastic disturbances in channel flows of viscoelastic fluids has recently 
been investigated using linear systems theory [12J. For the Oldroyd-B constitutive model, 
computations reported in Ref. |T2] demonstrated that streamwise-constant disturbances can 
be considerably amplified in the weakly inertial/strongly elastic regime. As in Newtonian 
fluids, this amplification is fundamentally nonmodal in nature: it cannot be described using 
the standard normal mode decomposition of classical hydrodynamic stability analysis O |6] . 
Rather, it arises due to an energy exchange involving the fluctuations in the streamwise/ wall- 
normal polymer stress and the wall-normal gradient of the streamwise velocity [T3] . 

The key parameters in viscoelastic fluids are: the viscosity ratio, (3 = rjs/ {rjs + rjp), where 
rjs and rjp are the solvent and polymer viscosities; the Reynolds number. Re = pUoL/ {rjs+rip), 
which represents the ratio of inertial to viscous forces; and the elasticity number, fi = We/ Re, 
which quantifies the strength of elastic forces relative to inertial forces. Here, We = XUo/L 
is the Weissenberg number, A is the fluid relaxation time, Uo is the largest base flow velocity, 
L is the channel half-width, and p is the fluid density. By modeling ambient disturbances as 
stochastic forcing to streamwise constant channel flows of Oldroyd-B fluids (with spanwise 
wavenumber kz), an explicit scaling of the variance (energy) amplification with the Reynolds 
number Re was developed in Ref. [13], E{kz]Re, [3, p) = f{kz;P,p)Re + g{kz;P,p)Re^, 
where / and g denote the i?e- independent functions. It is worth noting that E quantifies 
the ensemble-average energy density (associated with the velocity field) of the statistical 
steady-state [21 II]- In the same reference, numerical computations were used to conjecture 
that in flows with p ^ 1 the steady-state variance scales as 

Eikz;Re,P,p) ^ /(A;.;/3)i?e + g{kz; l3)pRe\ 

In this work, considering case of high /i, we apply singular perturbation techniques for 
systems with rapidly varying inputs ^14j to establish that the variance indeed depends affinely 
on p. 

The last expression should be compared to the expression for the variance amplification 
in Newtonian fluids [T^, E^{kz; Re) = fiy{kz)Re + giy{kz)Re^. At low Re the fc^-dependence 
of Ej^ is governed by fN^kz), and at high Re Ej^{kz) ~ Re^g^ikz) [15]. We note that the kz- 
dependence of both / and /a? is characterized by viscous dissipation, which clearly indicates 
that the behavior of Newtonian fluids with low Re is dominated by diffusion. On the other 
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hand, the (^-functions exhibit peaks at kz ~ 0(1); the values of kz where these peaks take 
place identify the spanwise length scales of the most energetic response to stochastic forcing 
in Newtonian fluids with high Re, and in Oldroyd-B fluids with high ^Re^. Our results thus 
show that-in viscoelastic fluids-even at arbitrarily low (but non-zero) Reynolds numbers 
the steady-state level of velocity variance may be dominated by the i?e^-term owing to the 
elastic amplification of disturbances. This reveals the subtle interplay between inertial and 
elastic forces in Oldroyd-B fluids with low Reynolds/high elasticity numbers [9| ITO]. 

We consider incompressible channel flows of Oldroyd-B fluids with e = <^ 1. The 
equations governing the dynamics (up to first order) of velocity (v = [ f i V2 fa ]^), pressure 
(p), and polymer stress tensor (r) fluctuations around base flow condition (v, r) are brought 
to a non-dimensional form by scaling length with L, velocity with Uo, polymer stresses with 
rjpUo/L, pressure with (77^ + rip)Uo/L, and time with A 

ev = -i?e(VvV + Vvv) - Vp + (1 - /3)V • r + pV\ + V^d, = V ■ v, 

r = Vv + (Vv)^ - T + ^(-VvT - VvT + r • Vv + r ■ Vv + (r ■ Vv)^ + (r ■ Vv)"^). 

Here, a dot signifies a partial derivative with respect to time t, V is the gradient, Vv = v- V, 
and Vi, V2, and ^3 are the velocity fluctuations in the streamwise (x), wall-normal (y), and 
spanwise {z) directions, respectively. If t denotes time normalized by the convective time 
scale L/Uo, then t = t/We. 

The linearized momentum equation is driven by the body force fluctuation vector d = 
[ di 3,2 d^Y ^ which is considered to be purely harmonic in the horizontal directions, and 
stochastic in the wall-normal direction and in time. This spatio-temporal body forcing will 
in turn yield velocity and polymer stress fluctuations of the same nature. Our goal is to 
study the steady-state level of variance maintained in v by assuming that d is a temporally 
stationary white Gaussian process with zero mean and unit variance. It is noteworthy that 
if d denotes the forcing in the standard non-dimensional momentum equation (with time 
scale L/Uo)-, then d(t) = \/We d{t). This body-force scaling is introduced to guarantee the 
same auto-correlation operators of d(t) and d(t) [T3]. 

We study the linearized model for streamwise-constant three-dimensional fluctuations, 
which means that the dynamics evolve in the (y, z)-plane, but flow fluctuations in three 
spatial directions are considered. This model is analyzed since the largest level of variance 
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FIG. 1: Block diagram of the linearized streamwise constant model. 

in stochastically forced channel flows of viscoelastic fluids is maintained by streamwise- 
constant fluctuations [12]. 

The linearized governing equations can be brought to an evolution form by removing 
the pressure from the equations and by expressing v in terms of the wall-normal veloc- 
ity/vorticity (f 2, 1-1^2) fluctuations [I2]. Furthermore, application of the Fourier transform in 
t and z allows for elimination of the polymer stresses from the evolution model, which results 
in an equivalent block diagram representation of the linearized system. The block diagram in 
FIG. [T] provides a systems-level view of the streamwise-constant model. The boxes represent 
different parts of the system and the circles denote summation of signals. Inputs into each 
box/circle are represented by lines with arrows directed toward the box/circle, and outputs 
of each box/circle are represented by lines with arrows leading away from the box/circle. 
The inputs specify the signals affecting subsystems, and the outputs designate the signals 
of interest or signals affecting other parts of the system |16j . 

All signals in FIG. [T] are functions of the wall-normal coordinate y, the spanwise wave- 
number kz, and the temporal frequency u, e.g. V2 = ^2(2/, k^^oj), with homogeneous Dirich- 
let boundary conditions on 002-, and homogeneous Cauchy (both Dirichlet and Neumann) 
boundary conditions on V2- The bold symbols in FIG. [l] denote the i?e-independent opera- 
tors. These operators act in the wall-normal direction and some of them are parameterized 
by kz (Fj and Gj, j = 1,2,3), while the others depend on kz, uj, /?, and e (Ji, J2, and 
Cp). The operators and Gj are given by {Fi = ik^, F2 = — A;^A"\ F3 = —ikzA^^dy}, 
{Gi = —{i/kz), G2 = I, G3 = {i/kz)dy}, and they describe the way the forcing enters into 
the evolution model, and the way the velocity fluctuations depend on the wall-normal ve- 
locity and vorticity. The operator Cp captures the coupling from V2 to 002, and it is deflned 
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as 

Cp = Cpi + Cp2/(e(l + M'), 

where Cpi = —ikzU'{y) denotes the vortex tihing term, and Cp2 = i^z(l — P)Cp2, Cp2 = 
U'{y)A. + 2U"{y)dy, denotes the term arising due to the work done by the polymer stresses 
on the flow. Finally, Ji and J2 govern the internal dynamics of the wall-normal vorticity and 
velocity fluctuations, respectively. These two operators are given by Jj = (1 + iuj)'Kj, where 
Kj = (e(icj)2l - {f3Tj - t\)\u - Tj)-i. Here, Ti = A and T2 = A^^A^ respectively, stand 
for the Squire and Orr-Sommerfeld operators in the streamwise-constant model of Newtonian 
fluids with Re = 1, 1 is the identity operator, A = dyy — A;^ is a Laplacian with homogeneous 
Dirichlet boundary conditions, A~^ is the inverse of the Laplacian, A^ = dyyyy — 2k'ldyy + k"^ 
with homogeneous Cauchy boundary conditions, i = v^— T, U{y) = y in Couette flow, 
U{y) = 1 — y"^ in Poiseuille flow, and U'{y) = dU{y)/dy. 

The velocity and forcing components are related by Vi = Hijdj, where Hjj denotes the 
ijth component of the frequency response operator H, v = Hd. From FIG. [T], it is clear 
that the H^- are determined by 

iiu{kz,uj;Re,p,e) = Hii(A;^, t^; /?, e), 

iiij{kz,uj;Re,P,e) = \^Re^ iiij{k^,uj; p, e), j = 2,3, 

B.ij{kz,uj;Re,p,e) = \^iiij{kz,uj; P,e), i,j = 2,3, 

iiii{kz,uj]Re,p,e) = 0, i = 2,3, 

where the Hjj represent the i?e-independent operators, {Hn = yeGiJiFi; Hij = 
Gi JiCpJ2Fj, j = 2,3; H^- = ^/eGi32Fj, i,j = 2,3}. The variance maintained in v 
is given by 

1 f°° 

E{kz) = — / trace (li{kz, u!)ii* {kz, uj)) du, 

where H* is the adjoint of operator H. (For notational brevity we have omitted the depen- 
dence on Re, (3, and e in the last expression.) Using the properties of the trace operator 
we have E{kz) = Yl'ij = i-^iji^z), where E^j is the variance maintained in Vi by stochasti- 
cally forcing the linearized model with dj. From the definitions of E and operators Hij, the 
following i?e-scaling of variance amplification is readily obtained 

E{K-Re,l3,e) = f{kz-l3,t)Re + g{k,- j3,e)Re\ 
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where / = /n + /22 + /23 + /32 + /sa, 9 = gi2 + gis, with fij{k^) = 
(l/27r) trace (Hjj(A;2, co')H*j(A;2, tu)) du; and similarly for gij. This expression for the 
steady-state variance of v is valid for all Re, /3, and e. Our objective is to establish how the 
i?e-independent functions / and g depend on e in viscoelastic flows with e = ^ 1. 

The Reynolds-number-independent contributions to the steady-state variance can be de- 
termined by recasting each Hjj in the so-called state-space form 

icij{y,kz,t) = Aij{k^)yiij{y,kz,t) + Bj{k^)dj{y, k^,t), 
Vi{y,k^,t) = Ci{k^):s.ij{y,kz,t), 

where Xj^ is a vector of state variables, and {dj,Vi) is the input-output pair for frequency 
response operator Hjj. The state-space models contain a set of first order evolution equations 
(for state variables) and a static-in-time equation describing the relation between the states 
and the outputs. It is a standard fact [3] that the variance of Vi sustained by dj is determined 
by trace (CjPjjC*) , where Pij denotes the steady-state auto-correlation operator of Xjj, 



which is found by solving the Lyapunov equation, Aj,Pjj -|- Pjj A*- 
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We now examine how functions fij depend on e. The i?e-independent operators {Hn 
Y^eGiJiFi; H.y = ^/eGi32Fj, i,j = 2,3} admit the following evolution equations 
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with {/c = 1 for z = 1; /c = 2 for i = 2,3}, homogeneous Dirichlet boundary conditions on 
Xn and Zu, and homogeneous Cauchy boundary conditions on Xij and Zij for i,j = 2,3. 
Eq. ([T]) is in the standard singularly perturbed form as the time-derivative of Zij is multiplied 
by a small positive parameter e puy . The approximate solutions for Xij and Zij can be found 
by performing a slow-fast decomposition of the system's dynamics. By setting e = in the 
Zjj-equation, we obtain Zij = — (^Xij + T'j^^Fjdj^ , which is not a valid approximation 
of Zij] this is because of the white noise component dj in Zij, which yields infinite variance 
of the difference between Zij and Zij irrespective of how small e is. Nevertheless, Zij can still 
be employed as an approximation of an input Zij to the Xjj-subsystem as the slow system 
filters out the white noise component in Zij, while an approximation of the fast subsystem 
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is given by e%j = pT^z^j + Fjdj, with Zij{t) = Zijj{t) - ^(t) + 0(e^/2) [H]. Thus, 

the slow component of Zij arises from a contribution of TkXij and not from a contribution 
of the white noise input dj, as would be common in singularly perturbed systems subject to 
slow inputs To summarize, the slow-fast decomposition of system ([T| is given by 
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where A,,-, = B,-, = -{1/P)T-'F„ A.,^j = (/3/e)T,, B,-^ = (l/e)F„ C,,, = 

— Gj, Cjj = Gj. It is easy to show that the steady-state auto-correlation operator 
of [ Xij^s Zij J takes the form 

(l/e)Z,, 

which yields fij = trace(CijZijC*j) + 0(e), with AijjZij + ZijA*jj = -BjjB*j. A more 
formal procedure [H] can be used to show /(fc^; /?, e) = /o(fc2; P) + XlfeLi fk{kz'-, P), e -C 1, 
where foih;f3) = trace(CijZijC*j), {i,j = 2,3; z = j = 1}. 

Next, we examine the e-dependence of terms responsible for the i?e^-scaling of the steady- 
state variance. From the definition of operators Ji, J2, and Cp it follows that Hij = 
(l/A/e)GiKi(e(l + iuYCpi + Cp2)K2Fj, j = 2,3. This yields an equivalent block diagram 
representation in FIG. [2] and leads to the following set of equations for Hij, j = 2, 3, 

= T2^p + i(3T2 -el)^ + F^dj, 
e0 = Ti0 + (/3Ti - el) + cp, v, = v^Gi0, 
<f = Cpi{ei> + 2e^ + e^) + Cp2^, 
with homogeneous Dirichlet boundary conditions on (p, and homogeneous Cauchy boundary 
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FIG. 2: Block diagram of Hij(A;^, /3, e), j = 2,3. 
where all operators are partitioned conformably with the elements of x and z, 
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Since output Vi in Eq. ^ does not depend on z, we conclude that the fast dynamics does 
not influence evolution of vi. The approximate dynamics of the slow subsystem are given by 
{xij-s = Aij^s^ij^s + Bj^Jj, Vl = A/T7eCi,^xij-s}, with Ci,^ = [ Gi ], Bj-^ = -A22^(0)B2 = 
[ - il//3)T,% f, and 



-A22'(0)A2l(0) 



-I 

-T;^ "'^Cp2 —I 



Thus, we have gij = (1/e) trace(Ci^s-PiiC* ^) + 0(1), where Pij denotes the auto- correlation 
operator of ^ij,s, Aij^sPij + A^^- = — Bj^sB* In fact, it can be shown that g{kz] P, e) = 
{^/^)ET=o^'9k{h-,P), e « 1, where go{h; f3) = ^ - = 2 trace(Ci,,Pi,C^,J. 

To recap, in streamwise-constant Poiseuille and Couette flows of Oldroyd-B fluids with 
/i ^ 1, the level of variance maintained in v is given by 



Re / 
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where fk and gk are functions independent of Re and fi. Therefore, for sufficiently large 
fi, E{kz] Re, P, /i) ~ Refo^kz] P) + fj,Re^go{kz; P), which implies that in elasticity-dominated 
flows, the terms responsible for the Re- and i?e^-scaling of the steady-state variance are, 
respectively, /i-independent and linearly dependent on /i. 

It is readily shown that fo{kz; P) = fo{kz)/P, where the base-flow-independent function 
fo{kz) is given by /o(fc^) = -0.5(trace (T^;^) +trace (T2 ^)). Moreover, go{kz; P) = go{kz){l- 
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FIG. 3: Plots of foikz) (a) and go{kz) (b); go{kz) in both Couette (solid curve) and Poiseuille 
(circles) flows is shown. 



(a) 



rV*f 1 . If ct 



-2 -10 12 



(b) 



10.2 



0.1 



-0.1 



-0.2 



FIG. 4: (color online) Streamwise velocity fluctuations vi(z, y) containing most variance in Couette 
(a) and Poiseuille (b) flows. 

/3)V/3, with go{kz) = (A;2/4) trace (T^iCp2T2'C;2T^i). Thus, for /x > 1, 



E{kz;Re,P,n) 



Re 



foikz) + fiRe\l-(3YUk 



which shows that the variance depends affinely on ^ and increases monotonically with a 
decrease in the ratio of the solvent viscosity to the total viscosity. 

The analytical expressions for trace (T~^), j = 1, 2, were derived in [15]; these can be used 
to obtain a formula for fo{kz 
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In Couette flow, 



fNikz), which is illustrated in FIG. 
the expression for ^o(^z) simplifies to —{k'^/A) trace (T^'^T^^), and an explicit fc^-dependence 
of Qo can be derived after some manipulation. The resulting formula for go{kz) is used to 



generate the plot in FIG. 3(b) from this plot we observe the non-monotonic character of 



go{kz), with peak values at kz ~ 2.07 (in Couette flow) and kz ~ 2.24 (in Poiseuille flow). 
Streamwise velocity flow structures that contain the most variance in respective flows with 
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these two spanwise wavenumbers are shown in FIG. |4j The most amphfied sets of fluctuations 
are given by high (hot colors) and low (cold colors) speed streaks. In Couette flow the streaks 
occupy the entire channel width, and in Poiseuille flow they are antisymmetric with respect 
to the channel's centerline. These flow structures have striking resemblance to the initial 
conditions responsible for the largest transient growth in Newtonian fluids [2]. Despite 
similarities, the fluctuations shown in FIG. |4] and in Ref. [2] arise due to fundamentally 
different physical mechanisms: in high i?e-flows of Newtonian fluids, vortex tilting is the 
main driving force for ampliflcation; in high /x/low -Re-flows of viscoelastic fluids, it is the 
energy exchange mechanism mentioned earlier [13]. 

In summary, a new mechanism for nonmodal ampliflcation of disturbances in elasticity- 
dominated flows is revealed using variance ampliflcation analysis of the linearized 
streamwise-constant model. This mechanism can serve as an initial stage of streak de- 
velopment which may-upon reaching a flnite amplitude-undergo secondary ampliflcation 
and/or instability and eventually lead to an 'elastic turbulence' P, [TU]. Elastic turbulence 
may flnd use in promoting mixing in microfluidic devices, where inertial effects are weak due 
to the small geometries [2]. 
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